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We define a mapping of the QCD Balitsky-Kovchegov equation in the diffusive approximation 
with noise and a generalized coupling allowing a common treatment of the fixed and running QCD 
couplings. It corresponds to the extension of the stochastic Fisher and Kolmogorov-Petrovsky- 
Piscounov equation to the radial wave propagation in a medium with negative-gradient absorption 
responsible for anomalous diffusion, non-integer dimension and damped noise fluctuations. We 
obtain its analytic traveling wave solutions with a new scaling curve and for running coupling a new 
scaling variable allowing to extend the range and validity of the geometric-scaling QCD prediction 
. . . beyond the previously known domain. 
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^ ' Geometric scaling is originally [1] an empirical scaling law satisfied by the data on deep-inelastic totally inclusive 
^— ^ , cross-sections. It may be translated as a property of the scattering amplitude 

T{L,Y)=T[r{L,Y)] (1) 
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where the scaling variable t(L, Y) = L — vY in the original formulation. The usual kinematic variables are Y, the total 
rapidity and L = log((3^/A^) with Q the photon virtuality and A the Quantum Chromodynamics (QCD) reference 
scale at one-loop level. 

r^' On the theoretical ground, there exists known encouraging features of QCD in the saturation regime leading to 
approximate geometric scaling, but the situation is not yet settled. The good news were, for instance, that in the 
large- TVc limit and in the mean-field approximation, we are led to consider the Balitsky-Kovchegov (BK) equation 
[2] which possesses the nice property [3] to be mapped, in the diffusive approximation, onto the Fisher-Kolmogorov- 
Petrovsky-Piscounov (F-KPP) equation [4, 5] . This equation which has been widely studied in statistical physics is 
known to admit traveling waves as asymptotic solutions which translate into geometric scaling in QCD [1]. In fact one 
may enlarge that result to a whole class of equations, including the one considering the full QCD kernel [6] beyond 
\^ the restriction to the first two derivatives. One may introduce the notion of "universality class" since the asymptotic 
C ■ properties of the solutions, including the approximate geometric scaling, do not depend either on initial conditions or 
^""^ on a precise form of the nonlinear part of the equation. 

^NJ ' However, many questions remain unsolved. On a theoretical ground, there are quite a few problems to ensure 
analytical scaling properties, even in the simpler case of BK equation with running coupling, see e.g. the discussion 
of Ref . [7] . Even more puzzling is the discovery that introducing a small cut-off [8] to the FKPP or BK equation or 
adding a small white noise [9], i.e. for a stochastic (s)FKPP or (s)BK equation, which after all reflects a physical 
constraint relying on the finiteness and randomness of the number of gluons, predicts a strong violation of geometric 
scaling not seen on data. On the other hand, the same stochasticity applied to the running coupling case seems from 
; I indirect evidence either to preserve geometric scaling [10] or to delay its violation to much higher energy [11]. 

I apart from the phenomenological problems, which we will discuss in a further publication, the aim of the present 

paper is on the general theoretical issue of geometric scaling from the BK nonlinear evolution equation, in particular 
for the running coupling case. Translating the problem in more mathematical terms, it seems that the "universality 
class" of the BK equation with runnning coupling has not yet been found or at least clarified. It is our goal to propose 
a solution to this problem, and explore its consequences which happen to give rise to an extended scaling property, 
as will become clear soon. In fact we will define the initial mathematical problem as the diffusive approximation of 
the (s)BK equation with generalized running coupling, extrapolating from the running or non-running case. 
One starts by considering the following one-parameter family 
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L" dyNiL, Y) = {A^dl + AidL + Ao} N{L, Y) - N^{L, Y) + ^kN[L,Y)v{L, Y) , (2) 

where the case n~l {resp. n — 0) corresponds to the logarithmically running (resp. non-running) coupling in the (s)BK 
equation within the diffusive approximation (1/i" has been divided on both sides of (2)). Universality properties. 
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independent from initial conditions and the form of the nonlinear damping, will provide the suitable identification of 
the universality class. In equation (2) 

x{-dL)-^A2dl + A,dL + Ao (3) 

represents for our purposes the expansion up to two derivatives of the linear kernel, either for the FKPP equation 
(where x{l) =7^ + 1) or for the diffusive approximation of the Balitsky, Fadin, Kuraev and Lipatov (BFKL) kernel 
[6] relevant for the BK equation. The strength of the noise term is specified by the phenomenological fudge factor k, 
and v{L,Y) is the standard white noise satisfying {v) = and {v{L,Y)v{L' ,Y')) = S(Y — Y')S(L — L') . 

Investigating the properties of the general nonlinear evolution equation (2), the paper is organized as follows. In 
section 2, we show that Eq.(2) can be exactly mapped onto the extension of the stochastic Fisher and Kolmogorov- 
Petrovsky-Piscounov equation for the radial wave propagation in a medium with absorptive negative gradient, non- 
integer dimension and damped noise fluctuations when n > 0. In section 3, we derive new scaling solutions for 
arbitrary n verifying the general recipe (1), defining for every value of n a scaling variable r(L, Y) and a new scaling 
curve T{t{L,Y)) valid in the forward front of the traveling wave in a region (in L) beyond but neighbouring the 
previously considered [3, 8] domain. By compatibility with those previous universal results, in section 4 we establish 
the universality forms, indexed by n, of the traveling wave solutions and focus on the physically interesting cases 
n = 0, 1. We check recovering known resiilts for n = 0, 1 and find a geometric scaling parameterization extended in 
a larger kinematic domain and for the running case a new scaling variable. Summary and outlook can be found in 
section 5. 



2. MAPPING 



The initial idea of our approach is to introduce a change of variables aiming at restoring the "standard" constant 
diffusion term, which determines the mathematical order of the (linear part of the) partial differential equation (2). 
In some sense one is guided by restoring the classical diffusion equation (or also "heat equation" ) dx = dx + • • • , 
where T is "time" and X is "space" . One then writes the Ansatz 

L = dL = ^ X'-^ dx , (4) 

with /3 to be determined for our purpose. 

After a straightforward transformation, Eq.(2) takes the form 

^ {^^^-^^ai + (a.I^X-^^^ + ^X^-^^^ dx + Ao^NiX, Y) - 

-nX,Y) + ^l^^^uiX,Y) (5) 

where the noise term has been conveniently renormalized such that {v{X,Y)jy(X' ,Y')) =S{Y — Y')S{X — X'). 

Now, asking for a constant diffusion coefficient can easily be achieved by equalizing the powers of the first two terms 
of (5) namely 

np = 2-2p^P= (6) 

leading to the equation 

{AoN{X,Y)-N^{X,Y)) + J N{X,Y) u{X,Y) . 



Introducing new variables 

Ao N 
T=^Y, andi7(X,T) = — 
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Eq. (7) gets mapped onto the "sFKPP- type" equation 
dU{X,T) 

with 



dj^UiX, T)+ ((l-/3)X-V7X-(i-«) dxU{X, T) + XX-^^^-'^^ (U-U^) + X-3(i-/3) U p{X, T) (7) 



{v{X, T)v{X\ T')) = 8{X-X')8{T-T') . (8) 

Note that, as weU-known, one may introduce an extra factor Vl — V in the noise term to ensure an absorptive 
boundary also at f7 = 1. However, the effects of the noise are really important only in the dilute tail of the wave, 
while this modification allows to recover exactly the Id sFKPP equation for /? = 1. 

A suggestive interpretation of Eq.(7) comes from the redefinitions (8). While the diffusion constant is by definition 
rescaled to 1, the successive coefficient functions acquire the following meaning: 

• {\ — ff)X~^dX is the radial Laplacian term corresponding to a medium of noninteger dimension 2 — /J. 

• 7X-(i-^) is, at /3 = 1, the "drift" term which is now decreasing at large X for /3 ^ 1 and thus no more absorbed 
by a redifinition of the "space" variable (see [3]). It may be interpreted as a kind of space-dependent rotation 
of the space-time frame. 



XX ''^ defines a dynamical birth-death coupling strength characterizing a medium equipped with a negative 
gradient; it has thus the important effect of a space-dependent "absorption" on the wave propagation. 

• /J, X^''^^^'') is the new dynamical noise strength, strongly decreasing when X gets large. 

A qualitative interpretation of the hierarchical behaviour of these coefficients as a function of X suggests a similar 
hierarchy of the expected effects on the solutions. Indeed, we will see that it leads to a /3-dependent anomalous 
effective diffusion. Moreover, the fluctuation strength (asssociated to the quadratic correlators of T)) leads to 
an expected damping of the fluctuations and thus of the stochasticity effects. In the rest of the paper we will thus 
neglect the stochasticity, postponing its precise study to further work. 

It is mathematically interesting to note that Eq.(7) possesses a well-defined limit n — > oo, i.e. /? — > 0, giving rise to 
the equation 

drUiX, T) = d%U{X, T) + X-^dxU{X, T) + \X-'^U{1 -U) + Vk X'^ v{X, T) , (9) 

which appears to be the radial sFKPP equation in a 2-dimcnsional medium characterized by a negative gradient 
corresponding to the birth-death factor \X^^. The noise term is also depressed by the factor kX~^ which defines the 
coupling strength of the noise correlator. 

Hence, in general, for nonzero /? the constant (by construction) diffusive coefficient is larger than the effective birth 
rate, which is itself larger than the effect of the noise correlators. This hierarchy can be interpreted as due to a 
"radial" sFKPP equation in a medium of dimensionality between one and two with negative gradient and decreased 
noise. Propagation of waves in such a "medium" should and will affect the properties of the traveling solutions which 
we shall now examine. 



3. SCALING SOLUTIONS 



Inspired by the fundamental traveling-wave property of the standard FKPP equation, our strategy will be to foc;us 
on some region of the forward front and thus look for a scaling variable t{L, Y) expressed in the new parameterization 
as X{X,T). It will lead to a solution of the form of a traveling wave U{X), based on an ansatz valid in a range of 
"space" and "time" which will be made precise later on. Let us recall that in all the following we shall concentrate 
our study on the deterministic equation without noise. 
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Scaling Solution: First Order Approximation 

Let us initially introduce an Ansatz 

U = Ae-^^-"^^"")' , (10) 

where the constants a, b, a and the "speed" v have to be determined. Our strategy is to first solving Eq.(7) at leading 
order when X X is large w.r.t. vTX~°' and in a second step going beyond this approximation. Using the convenient 
formulae 

drU = drX dxU dxU = dxX d^U dx^U = dx^X %t/ + [dxXf d\U. (11) 
we can write, within the same approximation X ~ X 

^drU = (abX^*-!) vX-"" 

^dxU = {-abX^''-^) {1 + vaTX"-'^) 

1-dlU = {a%^X^'~^ + ab{l-2b)X'''-^) {l + vaTX^-'f + 

+ {abX^''-^} {l-2h){2-2b)vTX^''-^ . (12) 

In a first step we solve the linear part of Eq.(7). As well-known [5], the forward front can be studied by restricting 
first to the solution of the linear equation since the nonlinear part has no direct effect (but important indirect selection 
effect) on the solution in this range. Selecting now leading orders when X X is large w.r.t. vTX~°' we find an 
iterative solution of (7) by ordering in decreasing powers of X. Within the approximation X ~ X and inserting the 
relations (12) in (7) one realizes that the leading order should be obtained with terms proportional to X'^^~'^. 
Collecting the leading terms from all contributions one gets 

I- drU : mbX^''-"-! 



U 



ux 



dxU : -ab{l-P)X 



2b-2 



lX^~^^ dxU : -^abX'^+^^-^ 

^d%U : aH^X^^-^ + ab{l-2b)X^^-^ 

^ XU : AX2^-2 . (13) 

One thus gets a matching of the dominant contribution in X'^l^~'^ by the following relations 

a=l-P; b=^ 
vab = a^b'^ - 706 + A (14) 
and thus finds two solutions for the parameter a, namely 

a± = ^ + 7 ± + 7)' - 4a) . (15) 

Hence, as expected from solutions of a second-order equation, the general solution can be (except for the case a+ = a_ 
see further on) expressed as 

U{X,T) - A+ e-'^-i^-^^^'-")^ - A_ e-«+(^-^^^-0^ . (16) 

Looking to the solutions (15) one sees that one has the condition v > 2^/\ in order to have real solutions, and thus 
the dominant speed is minimal when the condition 

Vc + ^ = 2\/A = ac/3 (17) 

is satisfied, where Vc (resp. a^) are the critical speed (resp. critical slope). We recognize among the conditions (17) 

the same critical speed as the one determined by the standard method [3, 8] when n = 0, 1 . Note however that the 
critical slope ac is /3-dependent contrary to what is found with the Brunet-Derrida type of Ansatz [3, 8] . 
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Scaling Solution: Second Order Approximation 

The critical conditions are realized as a special solution for the "degenerate case" a+ = a_ corresponding to the 
Ansatz (where we do not add the possible second pure exponential term, as will be justified in the following section) 

f/ = ^e-ac(x-.cTx'^-^)^ (^X-VcTX^-'y , (18) 
corresponding to a geometric scaling variable (restoring for a while the L and n dependence), 

X ^ X - UcTX^-i = X - VcTX-^ = - VcTL-^ . (19) 

Taking into account the prefactor in (18), the corresponding set of needed expressions (omitting from now on the 
subscripts v; ttc ^ a for simplicity) can be written using for simplicity the notation 



X 



{1-VcTX'^-^) ' . (20) 



One obtains 



OtU = U{al3/2X^-^- (5/2X^1^-^^ vX^-^w^-^ 

dxU = [/ (-a/3/2X'5-i + /3/2X'5/2-i) (i + v{l - (3)TXf^-^uf^-^) 

d%U = u(a(3/2X^-^~f3/2X^/^-^)' {l + v{l - (})TX^-^u^-^f + 

+ (7(a/3/2(l-/?)X'^-2-/3/2(l-/3/2)X'^/2-2) {l + v{l - l3)TX^-^uj'^-^f + 

+ U {a^l 2X^-^-13/ 2X1^ '^-^^ {l-(3){2-p)vTX^-^ujl^-^ . (21) 

Now, we find an iterative solution of (7) by ordering in decreasing powers of X. Assuming before further discussion 
to consider the kinematical region where 

uj = l+0 [TX^-^] ~ 1 , (22) 

and inserting the relations (21) in (7) order by order, and considering the leading terms from all contributions one 
gets 

^ OtU : va(3/2X^^-^ - t;/3/2X3/3/2-2 ^ . . . 

iy^axC/ : -a/J/2(l-/3)X'5-2 + ^/2(i_/3)x/3/2-2 + ... 

-tX"-^^ dxU : -7 (a/?/2X2/3-2 - /3/2X3/3/2-2^ ^ . . . 

^ dj,U : o2/3V4X2/3-2 _ 0^2/2x3/3/2-1 ^ ^^^/2(1 - /3) - /?V4) X'^"' - /3/2(l - (i/2)X'^/^-^ + ■■■ 

^ XU : AX2/3-2 . (23) 

Collecting the leading terms from all contributions to Eq.(7) one gets from the four first orders 

: a'^p'^/4-vap/2-jap/2 + X = 

t;/3/2 + 7/3/2 - a/3V2 = 
-a/3/2(l-/3) + a/3/2(l-/3) - = -(3^/4 

/3/2(l-/3) - /3/2(l-/3/2) = -13^4 . (24) 

Both first relations are satisfied by the critical condition (17) Vc + ^ = 2\/A = ac/3. Hence the consideration of 
the special solution (18) allowed for one subleading order vcryfying the equation. The nonzero remaining orders are 
Xl^~^,X20-'^. As we shall see further on, the subsubleading contributions, together with those coming beyond the 
approximation (22) will be dealt with the linearly T-dependent terms. Interestingly enough, these subsubleading 
powers will be related, in the known cases at least [3, 5], to the subleading terms of the universal expansion of the 
scaling variable. 



xi/3-2 
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4. EXTENDING GEOMETRIC SCALING 

In order to explore the kinematic range of validity of the scaling variable X{X, T) and look for corrections allowing 
to take into account sub-sublcading corrections induced by the set of equations (24) , let us recall the previously known 
results on geometric scaling coming for the studies [3] on nonrunning (n = 0) and running (n = 1) BK equations and 
make the appropriate generalization of these known results for generic n. 

One is led to distinguish 3 regions [5] (we will sec that we can add a 4th one) namely, in their order towards 
the moving tail of the wave, the wave interior region, where geometric scaling has been identified, the leading edge, 
where diffusive corrections appear, and the very forward region where the nonuniversal initial conditions prevail. 
Extrapolating for all n the known results for n = and 1, the wave interior region can be characterized by the 
kinematic relation L ~ (vcT)'^ and the leading edge region by L ^ (vcT)'^ + p{vcT)'^ where p ~ est. 

Focussing on the leading edge and extrapolating for all n the Ansatz proposed in refs. [3, 8] for the standard running 
and nonrunning FKPP cases, the generic solution reads 

Un{X,T) oc G„ f^ L-{v,T)^^ + c^{T) ^ g7{L-KT)^+c„(T)} ^ (25) 

where n = (0, 1) is the previously defined index for nonrunning and running cases, t„(L, T) = L — (vcT)'^ + Cn{T) 
is the original scaling variable; c;xtrapolated from the values obtained in [3], Cn{T) the corresponding subleading 
universal contributions to the scaling variable and 7 a universal critical slope value characterizing the scaling in the 
wave interior when z = Tn{L,T)/T'^ — > 0. Note the major constraint in this limit G„(z) ~ z which expresses 
the "absorptive" property [12] due to the nonlinear terms. This is the very same property which is realized in our 
Ansatz (18), Motivating the absence of a pure exponential term a priori allowed for a solution of a second-order linear 
equation. 

The function G„ has been determined for n = 0, 1 (see e.g. [3]), while their precise form in the general case, not 
relevant for our needs, should not be difficult to be obtained in a similar way. 

Within the typical leading edge region where L ~ (vcT)'^ — Cn{T)T + pnT^^ , where p„ is a smooth function 
around some constant value. We see that the expression (25) takes the scaling form 

Un{X, T) oc G„ ipn) e^P-^^ ~ — r„(L, T) Gn {pn) e^""(^'^) . (26) 

Pn 

Hence in the leading edge region, we can formulate the standard solution with essentially the same original scaling 
variable. 

Strikingly enough, a matching with our obtained solution (18) is obtained in an appropriate kinematic region where 

X ~ Kv,T ; L ^ {Kv,T) ^ , (27) 

where the constant if > 1 in the nonrunning case. This kinematic region is always near and forward to the identified 
leading edge regions characterized (see [3] for n = 0, 1) by i ~ VcT^^ + 0{T~^). For instance, L ~ (wc + K)T for 
the nonrunning case and L - v^T"^/^ for the running case. Let us call for convenience this region the extended scaling 
region [13]. Note also that one identifies 

ac = {n + 2)\/A ~ 7p„ . (28) 

Hence requiring a matching between the leading edge and extended scaling regions, and using the scaling variable 
X = X — VcTX"^ leads to an extended scaling domain. Interestingly enough the scaling variable gets modified, 
since 

X = X- VcTX-^ = - VcTL-^ . (29) 

As an application to the physical cases, in the nonrunning case, one realizes that the scaling variable is the original 
one, while the overall scaling region has been enhanced thanks to the modification of the scaling function. In the 

13 1 

running coupling case, one finds X — VcTX~s = — VcTL~^ , which shows the compatibility with the standard 
scaling variable t = L — (vcTY^'^, since ~ VcT near the wave front. 

We thus obtain a new scaling solution extended in the forward region matching at its lower X boundary with the 
leading edge. We thus infer that, for all value of /? = the variable X — VcTX^~^ = 1/7+^ —VcTL~^ is a promising 
candidate to ensure an extended geometric scaling. 
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In order to check the vahdity of the scahng in this domain, one has to ensure that the approximation (22) is vahd 
for X (X T. By inspection of the system of equations (21), one reaUzes that the correction is of order X^^~'^ << 
X2/3-2^2/3-3 at large X, ensuring that the two first relations of (24) are still valid, preserving the enhanced scaling 
property. Hence beyond the wave interior and the leading edge and before the very forward regions classically 
analyzed e.g. in Ref.[5], one finds a new extended scaling region for generalized FKPP equations and their traveling 
wave solutions. 

Concerning the subsubdominant orders in (24), they may receive contributions from those 3/3 — 3 terms coming 
from (22) (except for the nonrunning case /? = 1) and those in /? — 2 and ^ — 2. A simple way to take these terms into 
account is to add subleading terms in the T-dependence of the scaling variable taken into account by the function 
c(T) which appears in the leading edge formula (25). For instance, considering 

X = X-VoT + h\ogT + kT-^ + ■ ■ ■ for (3 = 1 ; X = X - VcT + hT^'l^ + kT'^ +■■■ for /? < 1 , (30) 

one generates the order [3 — 2 and f — 2. Other inputs will easily ensure the cancellation of the subdominant terms 
of order 3/3 — 3 present in the (3^1 cases. Note the interesting fact that the expressions (30) are exactly matching 
the orders of the well-known subdominant universal terms [3, 5]. In the case /3 ^ 1, new terms appear which deserve 
a specific study. We postpone this for further work. 



5. SUMMARY AND OUTLOOK 



Let us summarize our main results: 



Mapping to (s)FKPP universality classes. Wc introduce a one-paramctcr indexed by n family of nonlinear 
equations including and interpolating the Balitsky-Kovchegov equations with nonrunning and running QCD 
couplings as {logQ^}"". In the diffusive approximation, this family is exactly mapped to a (s)FKPP equation 
describing the propagation of radial traveling waves in an absorptive medium of dimensionality 2(n + l)/(n + 2) 
and coefficient functions depending on the "space" variable X = {log (3^}"+^. The different coefficients verify 
(except for n = 0) at large X the hierarchy 

X° (Diffusion) » X-^ {''Drift") » X'^ {Birth/Death rate) » X~(^ {Noise strength) (31) 

Geometric Scaling Solutions. Using an expansion at large X, we find a new geometric-scaling domain forward 
to the canonical traveling wave front such that L ^ {KvcT)'^ , where T oc Y the rapidity plays the role of 
"time" , Vc is the minimal speed (the constant K is arbitrary except that K > 1 when n = 0) with the scaling 
solution 



U(x e \ J (l^ -VcTL-^y , (32) 

where A is a constant depending on the BK equation parameters, see (8). 

• Extended Geometric Scaling. Matching with the universal solutions in the near-by less forward leading edge 
region, the scaling variable — VcTL~^ is shown to allow for an extension of the geometric-scaling region 

in a new kinematic domain. 

Specifying for the nonrunning case (n = 0), one has the same scaling variable as previously, but a new scaling 
form allowing to extend the geometric scaling region from the wave interior L ^ VcT, including the leading edge 
region with L (vcT) + poT^ up to an extended forward region where L ~ KvcT with K > 1. 

For the running and more general n 7^ cases, one finds starting from the wave interior region L ^ VcT^^ , 
scaling in the leading edge region L ~ (ticT) "+t -|- piT^^ up to an extended scaling region where L cx (I'cT) "+2 . 

The output of these properties is interesting in particular for the running coupling case {n = 1), which is physically 
meaningful for the BK equation in QCD. For instance, the parameter Oc = (n + 2)\/A related to an effective slope of 
the scaling curve is increasing by a factor | from nonrunning to the running case, while the leading edge form leads 
to a constant slope. It would be worth comparing this prediction with a simulation of the BK equation [14] where a 
similar effect has been noticed. 
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On the phcnomenological ground, which is not our subject in the present paper, there are quite a few issues which 
deserve to be studied. Last but not least, it would be worthwhile to confront the new extended scaling properties, in 
particular the new scaling variable X = X — v^TX^^ = — [vi.T)L~^ , with data using e.g. the method of Refs.[15]. 
Also, geometric scaling in other reactions, such as hard diffractive events [16], would be welcome. 

On the theoretical ground, some studies remain to be done in the new framework. For instance, the effect of a cut-off 
[8] could be analyzed using the general solution (16) with probably complex conjugate functions and match with the 
cut-off region. Also, the investigation of subleading universal terms is possible, thanks to appropriate modification 
of the time dependence of the scaling variable. It could be worthwhile to merge the new scaling approach with the 
studies on the parametric form of the traveling waves [17]. 

One would also like to go beyond the diffusive aproximation. At the present stage of our theoretical knowledge, 
the BK equation has been corrected for the full nest-leading log contributions [18]. Indeed, the matching with the 
standard universal formulation ensures that one should not expect dramatic effects, except for the neat difference 
between the fixed and running coupling cases [19], but this dedicated study deserves to be done in the future. In any 
case, it would be interesting to see whether and how it works. 
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